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ABSTRACT 

To  avoid  blind  speed  phenomena,  staggered  PRF  MTI  systems  are 
used  in  modern  radars  for  detection  of  high-speed  targets.  In  this 
report  a  procedure  is  presented  which  considers  both  the  effect  of 
filter  weights  and  the  interpulse  durations  simultaneously.  The 
design  goal  of  this  procedure  is  to  yield  a  filter  with  minimum 
variations  in  the  passband  and  maximum  attenuation  in  the  stopband. 
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DESIGN  OF  A  STAGGERED-PRF  MTI  FILTER 


INTRODUCTION 

An  MTI  (moving  target  indicator)  filter  commonly  known  as  a  comb 
filter  exhibits  a  periodic  property.  Its  frequency  response  folds  over 
at  the  radar  pulse  repetition  frequency  (PRF) .  Thus  a  target  having  a 
doppler  frequency  which  is  an  integer  multiple  of  the  radar  PRF  will  be 
treated  as  stationary  clutter  and  will  be  rejected  by  the  MTI  filter. 

This  blind  speed  is  a  function  of  the  radar  operating  frequency  and  the 
PRF  which  in  turn  determines  the  unambiguous  range.  For  a  modern  radar 
which  operates  in  the  microwave  range  (L-band  or  above)  and  which  is 
designed  to  have  a  large  unambiguous  range,  it  is  difficult  to  avoid  this 
blind  speed  problem,  in  particular  when  the  radar  is  required  to  detect 
targets  with  very  large  velocities.  One  of  the  approaches  to  alleviate 
this  problem  is  to  use  a  staggered-pulse  MTI  system  in  which  the  inter¬ 
pulse  durations  vary  from  pulse-to-pulse . 

The  design  of  such  filters  has  been  discussed  in  the  literature  (for 
example  in  references  T 1 , 2 , 3] ) .  However,  these  design  approaches  are 
concerned  only  with  the  filter  weights,  with  the  assumption  that  the 
interpulse  durations  are  fixed.  Furthermore,  most  other  procedures 
determine  a  set  of  weights  which  yield  a  good  clutter  rejection  in  the 
stopband  region.  The  ripples  and  variations  of  the  filter  response  in 


1 


2 


HSIAO  AND  KRETSCHMER,  JR. 


the  passband  are  generally  ignored.*  For  a  good  filter  performance,  the 
smoothness  of  the  response  in  the  passband  is  just  as  important  as  the 
attenuation  in  the  stopband.  That  is,  target  detectability  should  be 
uniform  in  the  passband.  Thus  the  existing  methods  solve  only  part  of 
the  problem.  A  desirable  solution  yields  a  filter  with  maximum  attenuation 
in  the  stopband  and  minimum  ripples  in  the  passband  region.  To  achieve 
this  kind  of  characteristic,  more  degrees  of  design  freedom  are  required. 

In  the  case  of  the  fixed-PRF  MTI  system,  these  are  provided  by  including 
feedback  loops  in  the  delay  line  network.  By  properly  adjusting  the 
locations  of  the  poles  and  zeros  of  the  transfer  function,  it  is  possible 
to  design  a  filter  having  specified  characteristics.  This  aspect  has  been 
adequately  treated  in  the  literature.  However,  feedback  loops  generally 
increase  the  filter  transient  response  time  which  is  undesirable  for  radar 
MTI  applications.  For  a  staggered-pulse  MTI  transversal  filter,  an  additional 
degree  of  design  freedom  is  provided  by  variation  of  the  interpulse  durations. 
Thus  at  least  in  theory,  one  should  be  able  to  use  this  additional  degree  of 
freedom  to  one's  advantage.  However,  achievement  of  an  optimal  solution 


*During  the  preparation  of  this  manuscript  a  paper  by  Jacomini,  (0.  J. 
Jacomini,  "Weighting  Factor  and  Transmission  Time  Optimization  in  Video 
MTI  Radar,"  IEEE  Trans  Aerospace  and  Electronic  Systems,  Vol  AES-8, 
pp  517-527,  July  1972)  considered  variation  of  the  interpulse  periods. 
However,  his  criterion  of  optimality  is  based  on  minimizing  a  linear  cost 
function  and  differs  from  the  criterion  presented  in  this  paper. 
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to  such  a  problem  is  extremely  difficult.  In  this  paper,  we  attempt 
to  develop  an  approach  for  solving  such  a  problem.  In  the  formulation 
of  the  optimization  procedure,  we  propose  a  condition  which  follows 
closely  the  criterion  of  a  conventional  comb  filter  in  which  the  attenuation 
in  the  stopband  is  maximum  while  variations  in  the  passband  are  minimum. 

A  solution  which  is  based  on  a  search  technique  is  presented  which  at  least 
yields  a  locally  optimal  solution  for  the  imposed  conditions. 

POWER  TRANSFER  FUNCTION  AND  AVERAGE  GAIN 

The  impulse  response  of  a  staggered  MTI  filter,  shown  in  Fig.  1,  can 
be  represented  as: 

N 

f(t)  =  Z  an  6 (t  -  Tn)  (1) 

n=0 


where 


n-1 

E  Ti,  T0  -  0 
i=0 


and  t ^  represents  the  interpulse  duration.  The  Fourier  transform  of  the 
above  expression  is: 


F(f)  -  I  ai  exp  (-j2jrf  Tt)  (2) 

where  f  is  the  doppler  frequency  of  the  radar  return.  The  power  spectrum 
density  of  this  transfer  function  (henceforth,  it  will  be  called  the  power 


*The  summation  indices  in  (2)  and  subsequent  equations  vary  from  0  to  N 
unless  stated  otherwise. 
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transfer  function) ,  is 


I F(f)  | 2  =  E  E  8i  a<  cos  2jtf(Ti  -  T<). 
i  j  J  J 

The  design  goal  is  then  to  find  a  set  of  weights  a^  and  a  set  of 
interpulse  durations  Ti  such  that  the  power  transfer  function  has  the 
desired  characteristics. 

Before  investigating  the  design  conditions  and  procedures,  it  will 
be  informative  to  first  discuss  some  properties  of  the  power  transfer 
function . 


a.  Since  we  are  only  interested  in  the  relative  amplitude  of  this 
function,  the  weights  can  be  normalized  by  any  constant  and 
the  characteristics  of  the  function  will  not  be  changed. 

b.  If  the  interpulse  durations  are  constrained  such  that 

fr  Ti  =  ii 

where  ii  are  integers  for  all  i  then  fr  is  the  first  blind 
doppler  frequency.  This  is  equivalent  to  saying  that  the 
frequency  response  of  the  filter  is  a  periodic  function  having 
a  period  of  fr,  or 

| F(f ) | S  *  | F(f  +  nfr)|3. 

This  becomes  evident  if  one  replaces  f  in  (3)  by  i  +  nfr. 

For  the  case  that  the  interpulse  durations  are  fixed  having 
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Therefore,  in  this  case,  the  first  blind  doppler  frequency 
occurs  at  the  PRF  which  is  the  reciprocal  of  the  interpulse 
duration  t. 

c.  The  power  transfer  function  is  synmetric  within  each  period, 
that  is 

| F (f r  -  f)|S  =  E  Z  ai  aj  cos  2jt  (fr  "  f)  (Ti  "  Tj)  =  lF(f)|2- 

d.  The  average  value  of  this  power  transfer  function  is 

fr 

|F(f)|*  -  -r-  J  E  E  aj  a.  cos  2jt  (Tj.  -  Tj)  df  =  E  . 

rr  o  i  j  f 

e.  If  one  assumes  that  white  noise  having  a  unit  spectral  density 
function 

n0  -  1 

is  applied  to  this  filter,  the  ensemble  average  output  noise 
power  is  equal  to  ^  ai*.  Thus  if  the  filter  weights  are 

normalized  by  a  factor  J C  the  average  output  is  just 
unity.  Since  both  the  average  gain  of  the  power  transfer 

g 

function  and  the  noise  power  gain  are  equal  to  E  a1  ,  the 

i 

average  S/N  is  not  degraded.  Thus,  one  may  define  this  as 
an  average  gain  of  the  power  transfer  function. 

CLUTTER  OUTPUT 

Assume  that  the  power  spectrum  density  function  of  the  clutter 
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return  is  G(f),  then  the  clutter  output  power  is 

CO 

c(Ti>ai)  =  J  S  E  ai  aj  G(f)  cos  2jt  f  (Tt  -  Tj)  df. 


Changing  the  order  of  summations  and  integration,  one  finds 


C(Ti>a^)  =  E  E  a^  aj  p^j 

where 

CO 

P±j  =  J  G(f )  cos  2ir  f  (Ti  -  Tj)  df  , 

-GO 


(10) 


(11) 


which  is  seen  to  be  the  clutter  covariance  function  at  time  Tj.  -  Tj . 

One  of  the  most  desirable  optimization  conditions  is  to  minimize  this 
clutter  power  output.  However,  if  one  attempts  to  do  so,  it  will  lead 
to  a  trivial  solution  that  aj  *  0.  To  avoid  this  situation,  one  may 
require  that  the  average  gain  of  this  clutter  output  be  some  fraction  of 
the  average  gain  of  the  power  transfer  function  such  that 

?  ?  ®i  aj  Pi j  =  X  E  a^  (12) 

1  j  1 

where  X  represents  the  clutter  suppression  factor.  In  order  to  find  a 
set  of  aj  with  a  given  set  of  Tj  such  that  X  is  minimum,  the  following 
conditions  must  be  satisfied 

E  aj  p  -  X  at  »  0  (13) 

j  ij 

i  =  0,  ....  N 
j  *  0,  . . . ,N. 

This  is  a  typical  eigen-value  problem.  There  are  (N  +  1)  such  a^  sets 
which  satisfy  the  above  equation.  Each  of  these  sets  comprise  the 
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eigen-vector  associated  with  an  eigen-value  X.  Thus  the  eigen-vector 
associated  with  the  minimum  X  is  the  desired  optimal  solution.  In 
references  [l]  and  [3],  the  average  S/C  gain  of  the  MTI  filter  is 
optimized.  This  is  referred  to  as  the  reference  gain  and  is  numerically 
equal  to  the  reciprocal  of  X.  Alternatively,  this  procedure  is  equivalent 
to  maximizing  the  clutter  attenuation  of  a  filter  normalized  i.o  the  average 
power  gain. 

As  an  example,  assume  that  G(f)  has  a  uniform  distribution  function 
with  a  limited  bandwidth  such  that 


G(f)  =  - - - i  s  f  ^  ^u 


fu  *  fi 


0,  otherwise 


u 

- -  cos  2jt  f(Ti  -  Ti)  df 

-  f#  if 


ru  '  H 

sin  2rt  fu(Tt  -  Tj)  -  sin  2n  tj  (Tt  -  Tj) 

2*  (fu  -  f/)  (Ti  -  Tj) 

If  f£  »  0,  this  then  becomes 

sin  2t  fn(Tt  ■  Tj) 
lJ  2k  fu(Ii  -  lj) 

If  the  G(f )  is  Gaussian  and  has  a  mean  f0  and  variance  oc  »  then 
-  cos  2jt  f0  (Ti  -  Tj)  exp  [  -  2jtS  (Ti  -  Tj)3  •  oc3  ]  • 

In  the  case  f0  =  0,  this  becomes  the  familiar  case  treated  by 


Emerson  [l]  and  Murakami  [3]* 
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AVERAGE  GAIN  IN  THE  PASSBAND  REGION 

Next,  consider  the  stopband  and  the  passband  regions  of  the  power 
transfer  function.  Let  the  stopband  region  extend  from  f  =  0  to  fu  and 
the  passband  region  extend  from  f  =  fu  to  fr  -  fu.  Then 


I 


E  E  8j  cos  2jr  f  (Ti  -  Tj) 


df 


U 

2  J  E  E  fij  cos  2rt  f  (Ti  -  Tj)  df  + 


fr/2 


2  J  E  E  aj  8j  cos  2*  f  (Tt  -  Tj)  df. 


(18) 


Formulation  of  (18)  uses  the  fact  that  the  filter  power  transfer  function 
is  symmetric  about  fr/2. 

The  last  term  in  the  above  equation  expresses  the  output  of  a  signal 
which  has  a  uniform  spectrum  in  the  passband.  This  merely  states  that 
the  doppler  shift  of  the  return  target  signal  is  not  known  a  priori. 
Furthermore,  the  first  term  in  (18)  is  equal  to  E  at*  from  (8),  irrespective 
of  the  value  of  ai«  Under  this  condition,  the  average  output  in  the  pass- 
band  region  is  determined  from  (18)  as 

f  /  2 


- i -  |  T  Z  Z  a.  a.  cos  2*  f  (T,  -  T.)  df 

fr/2  -  fu  f  i  J  1  j  1  j 


fr/2 


fr/2  -  fu  1 


f 

.  E  n8  -  —1 -  J  ES«i»j  cos  2*  f  (^  -  Tj)  df.  (19) 


fr/2  -  fu  o  1  j 
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By  regarding  the  a^  to  be  derived  for  a  uniform  spectral  density  function 
of  unit  amplitude  from  0  to  fu  and  by  noticing  the  relations  of  (16)  and 
(13) ,  one  finds  that 


V2 

* 

,  EE  aA  cos  2jr  f  (T  -  Tj)  df 
^u  ^ 


-  [  1  + 


2fu  (1  -  X) 
fr  -  2fu 


]  S  *iS 


(20) 


The  left  side  of  this  equation  represents  the  average  response  of  this 
filter  in  the  passband  region.  In  the  case  fu  <  <  fr,  which  usually  is 
the  case  when  the  clutter  spectrum  is  concentrated  in  the  region  close 
to  f  =  0,  the  left  side  of  the  above  equation  can  be  approximated  by  the 
average  gain  ^  a^s .  Although  (20)  was  derived  for  a  uniform  spectral 
density  over  the  stopband  region,  the  effect  of  any  other  spectral 
density  is  to  modify  X.  For  any  practical  filter,  X  <  <  1,  so  that  the 
above  relation  can  be  generalized  to  include  all  other  cases. 

The  mean-squared  deviation  from  this  average  value  in  the  passband 
is  then 

V2 

era  « - 1 -  |  [  E  E  ai  a,  cos  2*  (Ti  -  Tj)]a  df.  (21) 

ir/2  -  fu  c  i*r j 

u 

In  order  to  form  a  filter  that  has  a  response  curve  with  minimum 
variations  in  the  passband  region,  this  error  function  must  be  a  minimum. 
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In  order  for  this  to  be  achieved  by  choice  of  a  set  of  T^,  it  is  required 
that: 


Then,  at  least  in  theory,  with  the  condition  of  (22)  and  (13),  one 
should  be  able  to  solve  for  a  set  of  and  a^  which  yields  a  filter 
design  having  the  desired  characteristics;  i.e.,  minimum  ripple  in  the 
passband  and  maximum  attenuation  in  the  stopband. 

Unfortunately,  the  condition  of  (22)  is  not  readily  solvable  in 
closed  form  by  known  methods.  In  the  following  section  we  propose  an 
approximate  means  to  solve  this  problem. 

A  SEARCH  SOLUTION 

In  view  of  the  difficulty  in  seeking  a  solution  for  the  above  problem, 
we  suggest  here  an  alternative  means  to  seek  an  approximate  solution. 

This  method  is  based  on  a  search  alogrithm. 

The  procedure  of  this  approach  is  shown  in  Fig.  2.  The  index  i 
represents  the  i1-*1  component  of  while  index  k  represents  the  success 
of  each  trial.  The  aj's  represent  the  filter  weights.  Index  M  represents 
the  number  of  iterations.  First  an  initial  set  of  interpulse  durations 
is  assumed.  With  the  assumed  interpulse  durations  and  by  use  of  the 
relation  of  (14),  a  set  of  optimal  weights  with  an  associated  eigen-value, 
X,  are  found.  If  this  X  is  less  than  the  specified  clutter  attenuation 
(S.C.A.),  these  weights  are  inserted  into  (21)  to  compute  the  value  of 
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the  error  function.  Both  the  eigen-value  and  the  error  function 
are  recorded.  The  next  step  is  then  to  either  increase  or  decrease  the 
interpulse  durations  by  an  amount  At^.  This  is  done  one  interpulse- 
duration  at  a  time.  Again  the  eigen-value  X  is  computed  which  is  then 
compared  with  the  specified  clutter  attenuation  factor  (S.C.A.)  and  the 
previous  recorded  X.  If  it  is  improved,  then  the  error  function  e  is 
computed  and  the  result  is  compared  with  the  previous  recorded  e.  If 
e  is  improved  then  the  new  value  is  retained  and  another  of  the 
interpulse  spacings  will  be  varied.  This  procedure  is  repeated  until 
no  further  improvement  in  e  is  obtained  and  the  function  has  apparently 
reached  a  local  minimum.  The  procedure  can  be  either  stopped,  or  a 
new  set  of  initial  values  of  can  be  assumed  and  the  same  procedure 
repeated  again.  It  is  assumed  that  an  adequate  number  of  pulses  are 
available  to  achieve  the  required  clutter  attenuation. 

Equation  (4),  repeated  below,  sets  the  boundary  conditions  on  the 
interpulse  duration 

fr  T i  = 

where  is  the  interpulse  duration  and  l ^  is  an  integer.  Hence,  the 
variations  of  t ^  are  limited  to  discrete  values  at  steps  which  are 
the  reciprocal  of  fr.  In  practice,  in  order  to  obtain  sufficient  energy 
on  a  target,  the  interpulse  time  variations  are  bounded  and  therefore 
only  a  few  values  of  need  be  considered  in  using  the  search  program. 
Moreover,  the  minimum  interpulse  time  is  dictated  in  practice  by  the 
maximum  range  of  the  clutter.  The  constraint  which  requires  that  T ^ 
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be  discrete  does  not  necessarily  satisfy  the  condition  of  (22).  By 
relying  on  the  search  method  discussed  above,  we  may  find  a  solution 
which  is  not  necessarily  absolutely  optimal;  however  the  solution  is 
always  physically  realizable.  Furthermore,  this  procedure  can  be 
repeated,  each  time  with  a  set  of  different  initial  values,  until 
satisfactory  results  are  achieved.  In  the  following  section  we  present 
several  examples. 

NUMERICAL  EXAMPLES 

All  examples  are  for  7-pulse,  staggered-PRF,  MTI  cancellers.  To 
reduce  local  minimum  problems,  the  search  routine  has  been  repeated  about 
ten  times.  Each  time  a  different  set  of  initial  interpulse  periods  are 
used.  The  best  results  are  then  plotted  as  a  frequency  response  curve 
shown  in  these  illustrations.  The  initial  T £  values  are  generated  by  a 
random  number  generator  as  discussed  in  the  previous  section.  The  whole 
search  routine  takes  about  two  minutes  on  a  CDC  3800  machine. 

Figure  3  shows  a  design  with  a  stopband  extending  from  f  *  0  to 
f  =  .2  of  fj  where  f<j  is  a  normalizing  factor  corresponding  to  the 
reciprocal  of  the  smallest  interpulse  period.  This  is  equivalent  to  a 
filter  notch  of  approximately  120  knots  for  an  L-band  radar  at  an  unambiguous 
range  of  80  nautical  miles.  The  firjt  normalized  blind  speed  occurs  at 
f  *  10.  The  eigen-value  which  represents  the  clutter  suppression  factor 
is  about  0.00004  and  the  average  rras  error  in  the  passband  is  0.670768. 

The  interpulse  durations  shown  in  that  figure  are  normalized  with  respect 
to  the  smallest  t^,  while  the  filter  weights  are  normalized  to  the  square 


•w : 


—  • 


NRL  REPORT  7545  13 

root  of  the  average  gain.  At  zero  dB  level,  the  gain 
of  the  power  transfer  function  is  equal  to  the  average  gain.  In  the 
second  example  shown  in  Fig.  4,  the  stopband  is  set  at  a  range  from 
f  =  0.5  to  f  =  0.657.  This  curve  clearly  demonstrates  the  periodic 
property  of  this  type  of  filter.  In  these  examples,  it  is  assumed  that 
the  input  clutter  power  density  function  is  uniform,  that  is 

G(f )  =1  ffSf^fu 

=  0  otherwise. 

However  this  same  procedure  can  be  applied  to  a  clutter  power  density 
function  with  an  arbitrary  distribution. 

SUMMARY 

The  first  order  effects  of  varying  the  interpulse  periods  and  weights 
in  an  MTI  transversal  filter  are  to  alter  the  filter  response  in  the 
passband  and  rejection  regions  respectively.  The  purpose  of  this  paper 
is  to  show  how  an  acceptable  filter  design  which  is  compatible  with  a 
radar  waveform  requirement  may  be  achieved  by  systematically  varying  the 
interpulse  periods  using  a  search  algorithm  to  minimize  the  response 
variation  in  the  passband  or  visibility  region  of  the  filter  while 
simultaneously  selecting  a  set  of  weights  which  maximize  the  clutter 


attenuation. 
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Fig.  2  —  Flow  chart  of  aearch  program 
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Fig.  3  —  Seven-pulse  staggered  PRF  MT1  filter  frequency  respunse  curve  (fk  -  0,  fu  -  0.2; 
N  -  7;  X  -  0.000040;  C,  -  0.670768;  r,  -  1.1,  1.1,  1.0,  1.2,  1.4,  1.1;  a,  -  0.080892, 
-0.301091,  0.584633,  -0.617309,  0.379645,  -0.179466,  0.059962). 
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Fig.  4  —  Seven-pulse  staggered  PRF  MTI  filter  frequency  response  curve  (fv  «  0.5,  fu  =  0.657; 
N  -  7,  X  -  0.00004;  e,  »  0.624926;  T,  -  1.3,  1.2,  1.0,  1.1,  1.3,  1.1;  a,  •  0.075466.  0.120090, 
0.556937,  0.716912,  0.342963,  0.177940,  0.080346). 
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